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Abstract. The electronic and vibrational contributions to the static dipole 
polarizability of C70 fullerene are determined using the finite-field method within 
the density functional formalism. Large polarized Gaussian basis sets augmented 
with diffuse functions are used and the exchange-correlation effects are described 
within the Perdew-Burke-Ernzerhof generalized gradient approximation (PBE-GGA). 
The calculated polarizability (a) of C70 is 103 A'' , in excellent agreement with 
the experimental value of 102 A"^ and is completely determined by the electronic 
part, vibrational contribution being negligible. The ratio a(C7o)/a(C6o) is 1.26. 
The comparison of polarizability calculated with only local terms (LDA) in the PBE 
functional to that obtained with PBE-GGA shows that LDA is sufficient to determine 
the static dipole polarizability of C70. 
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C70 is perhaps tlie most studied carbon fullerene after the Cgo fullerene. Several studies 
have addressed electric response properties of C70 and its more famous cousin Ceo 
[Il[a[3lia[5l[6l[71[8l[9l[T0lIIliaE3inil5l[T6l[I3 Static dipole polarizabihty ao ofCyo, 
which is the subject of present work, was recently measured in gas phase by Compagnon 
and coworkers[3]. Using the molecular beam deflection technique, they reported the 
mean static dipole polarizability to be 102 with an error bar of ±14 A^. Theoretically, 
polarizability of C70 has been subject of several investigations [U El El [9], [151 IE]- Only 
two of these studies are, however, at the ab initio level of theory. The first study is 
by Jonsson and coworkers[5], who using the self-consistent-field (SCF), and the multi- 
configuration self-consistent field (MCSCF) theories in combination with the 6-31++G 
basis set reported polarizability of C70 to be 89.8 A^. The second study is due to van 
Fassen et al. [16] who used time dependent density functional theory (TDDFT) and 
three different basis sets to calculate ao of C70 p^. Using the largest basis (triple 
zeta with additional field induced polarization function TZVP+) they found ao to be 
104.8 A'^. These authors also computed cto using the current-dependent Vignale-Kohn 
functional and concluded that polarizabilities calculated using the current- dependent 
functional gives good agreement with experimental values for Ceo and C70 . The 
geometric structure of C70 was not optimized in both these works. 

The present article complements these earlier studies and reports the static 
dipole polarizability of C70 calculated within density functional formalism, using large 
polarized Gaussian basis sets augmented by diffuse functions. The calculations are 
performed within the generalized gradient approximation using Perdew-Burke-Ernzerhof 
parametrization[18]. Unlike in previous works, we first determine the equilibrium 
structure of the C70 fullerene at the same level of theory. We then compute the 
static dipole polarizability. We investigate the vibrational a^b contribution to the 
polarizability cto as well as the electronic a^i contribution. The former, which has not 
yet been computed, is computationally significantly more expensive than the latter as 
it requires multiple optimization of C70 geometry in the presence of electric field or 
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Figure 1. The optimized geometry of C70 fullerene. The five inequivalent atoms 
in the D^ii structure are denoted by lower case letters. The bond distance between 
labeled atoms are given in Table I. 

requires the calculation of full vibrational spectrum. 
2. Computational details and Results 

The geometry of C70 was fully optimized using the limited memory Broyden-Fletcher- 
Goldfarb-Shanno algorithm [19] using the NRLMOL suit of codes developed by Pederson 
and coworkers [20], [211 122]. The code employs an efficient scheme for numerical 
quadrature for the exchange-correlation integrals [23]. The D^h symmetry[24j of 
C70 fullerene was exploited to reduce the computational expenses during structure 
optimization. The molecular orbitals in NRLMOL are expressed as a linear combination 
of Gaussian orbitals. The Gaussian basis set for C consists of 5 s—, 3 p—, and 1 d— 
type Gaussians each contracted from 12 primitive functions. The total number of basis 
functions used for the structure optimization is 2450 and that used for polarizability 
calculations is 2870. The exponents in the basis set are optimized iteratively by 
performing a self-consistent calculation on isolated atoms[22]. More details about basis 
set and its construction can be found in Ref. [22112S112S]- The fully optimized structure 
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of C70 is shown in Fig. 1. There are five inequivalent atoms in C70. The optimized 
positions of these inequivalent atoms in atomic units are a ( -2.3127, 0.0000, -7.4352 
), b ( -4.5356, 0.0000, 6.0395 ), c ( 5.5080, 1.3063, 4.5799 ), d ( 5.9336, 2.6454, 2.2646 
), e ( 6.5245, 1.3720, 0.0000 ). The C70 fullerene geometry can be generated using the 
xyz coordinates of inequivalent atoms and symmetry operations of point group D^h, 
where the five-fold highest symmetry axis is the 2— axis. In Table I, the calculated 
bond distances are compared with some experimental and theoretical values reported in 
literature. Agreement with experimental bond distances is quite good. The equatorial 
bond distance of 1.45 A is smaller than the 1.538 A measured in the gaseous electron 
diffraction experiment. However, the smaller values comparable to the prediction of the 
present calculation, have been reported in other experiments. The distance between the 
polar pentagons is 7.869 A in good agreement with experimental value of 7.906 A from 
the most recent experimental measurement |24j . 

The elements of the static dipole polarizability tensor are 

dF.dF/ 

Here, E is the total molecular energy and Fj is the i*^ component of the electric field. 
Eq. ([T]) is the coefficient of the second term in Taylor's expansion of total energy E in 
the presence of field: 

Alternatively, the polarizability tensor elements could also be obtained from the induced 
dipole moments. A number of methods have been formulated to obtain aij and several 
review articles describing the details and applications of these methods exist |27[ 128] . 
In this work we compute numerically using finite-difference formula. This requires 
calculation of total energy of the molecule for various field values. The electronic 
contribution to the polarizability agi is then obtained using a suitable approximation 
to compute the second derivative in Eq. [H Alternatively, least-square fitting technique 
can also be used to extract the polarizability j29| [30l [TO] . This way of calculating agi is 
known as the finite-field (FF) method. The implementation of FF method is fairly 
straight-forward. The Hamiltonian of the system is augmented by the term —F ■ r that 
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represent interaction of an electron with applied electric field F. The self-consistent 
solution is performed to obtain total energy for a given value of electric field. The 
self-consistent process takes into account the field-induced polarization or the so-called 
screening effects. The drawback is that addition of — F • r to the Hamiltonian, in 
general, breaks the full point-group symmetry. Thus computational advantage of the 
full point group symmetry is lost. In some cases, it may be still possible to use lower 
symmetry for applied fields along a specific direction. The present calculations of aei of 
C70 did not use any symmetry. The electric field step was chosen by ensuring that the 
calculated polarizability is accurate and uncontaminated by higher polarizability. The 
aei obtained from induced dipole moments agrees with that obtained from energy within 
1%. The basis set used during structure optimization was augmented by a d— type 
diffuse function with 0.0772097 exponent. This methodology has been found to provide 
a good description of aei of Ceo and several molecules and clusters [ini [HI EH [251 [26] . 
For more details of methodology and application we refer reader to a recent review by 
Pederson and Baruah[2B]. The computational scheme used in this work has also been 
applied to Cgo fullerene. In fact, one of the earliest calculation of polarizability of 
Ceo by Pederson and Quong used the same methodology [TT] . As mentioned earlier, 
the nuclear positions were assumed to be frozen during the calculation of a^i ■ The 
relaxation of nuclear positions in the presence of applied electric field also contribute to 
the polarizability. This contribution is often called vibrational polarizability amt, and 
usually is the second largest contribution of the polarizability. In some cases, particularly 
for the system with ionic or hydrogen bonding it could be even larger than aei ■ Here, we 
compute ayih within double-harmonic approximation. A full account of the formulation 
of calculation of a^b used in this work can be found in Ref. 

The mean or average polarizability aei , is one-third of the trace of the polarizability 
matrix: 

o^m = ^{axx + ctyy + a^z)- (3) 

The calculated values of the polarizability components are a^x = c^yy = 99A^ and 
azz = lllA'^ , where z— axis is along the five- fold axis. The mean polarizability 
aei is 103 A'^ . We have also calculated unscreened polarizability using the sum-over- 
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states expression with excitation energies approximated by eigenvalue differences. As 
in case of Ceo fuUerene [TOl [H], tlie unscreened polarizability turns out to be 330 
, roughly three times larger than the screened polarizability. In Table 11 calculated 
aei is compared with earlier published theoretical and experimental values. Recent 
measurements of gas-phase polarizability of C70 indicate aei to be 102 A^ with a 
rather larger error bar of about 14 A^ . The predicted (PBE/NRLMOL) value of 
103 A^ agrees well with the experimental value. As can be seen from the Table 11, 
all ab-initio values fall within the error bar. The PBE/NRLMOL Oei agree well with 
time dependent density functional theory calculation with statistically averaged orbital 
potential (TDDFT/SAOP) polarizability [16]. But, it is larger than the predictions 
by Hartree-Fock (HF) theory (HF/6-3H-G)[32j and time dependent current density 
functional theory (TDCDFT/VK)[16jmodels. The differences in the structure of C70 
used in those calculations and that optimized in this work is one source of discrepancy. 
In particular, the higher order polarizabilities are often sensitive to the molecular 
structure [33J. The other possible cause of differences is in the modeling of many body 
effects. In our recent study that compared static dipole polarizabilities of 142 small 
molecules in the HF and PBE models, the HF polarizabilities were found to be smaller 
than their PBE counterparts |34j . The larger value of cVe/in the PBE/NRLMOL than in 
the HF model is consistent with this observation. The aei obtained within the time- 
dependent current density functional model (TDCFT/VK) falls intermediate between 
the HF/6-31+G and PBE/NRLMOL predictions. We repeated calculations with only 
local terms in the PBE functional, to estimate the effect of gradient correction to the 
exchange-correlation functional. This approximation gives aei to be 100 , indicating 
that local approximation is sufficient to determine the dipole polarizability of C70. 
The vibrational contribution to the polarizability tensor within the double harmonic 
approximation [25] is given as amb,i,j = ■^'j^^jCj^^Zj^. Here, uj^ is the frequency of the 
/ith vibrational mode, Zj^^ is the effective charge tensor (See Ref. [25] for details). The 
vibrational contribution along the five- fold axis is a'^f is 0.43 A'^ while that along 
the transverse axis is 0.74 A'^ . Thus, in comparison with electronic polarizability, 
vibrational contribution to the dipole polarizability of C70 is negligible and the static 
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Table 1. The comparison of selected bond distances calculated in this work 
(PBE/NRLMOL) with those reported in literature. GED is the gaseuous electron 
diffraction measurements [24 ; Solid-State Electron diffraction (SED)[35]; Neutron 
difftraction (ND) mcsurmcnts [36]; X-ray diffraction [S^; Hartree-Fock/Dobule Zeta 
Basis [38]; BP86/TZP[3|; PBE/NRLMOL (present). 



Bond distance 


GED 


SED 


ND 


X-ray 


SCF/DZP 


BP86/TZP 


PBE/NRLMOL 


C1-C2 


1.461 


1.464 


1.460 


1.458 


1.451 


1.454 


1.439 


C1-C6 


1.388 


1.37 


1.382 


1.380 


1.375 


1.401 


1.389 


C6-C7 


1.453 


1.47 


1.449 


1.459 


1.446 


1.450 


1.436 


C10-C12 


1.386 


1.37 


1.396 


1.370 


1.361 


1.395 


1.383 


C7-C8 


1.468 


1.46 


1.464 


1.460 


1.457 


1.449 


1.433 


C8-C9 


1.425 


1.47 


1.420 


1.430 


1.415 


1.441 


1.426 


C9-C31 


1.405 


1.39 


1.415 


1.407 


1.407 


1.424 


1.410 


C31-C32 


1.538 


1.41 


1.477 


1.476 


1.475 


1.471 


1.452 



dipole polarizability of C70 is completely determined by the electronic polarizability. 
This result in consistent with the finding of Pederson and coworkers |25j for the Cgo and 
is relevant for the dielectric response of carbon nanotubes. 

To conclude, static dipole polarizability of C70 is calculated within density 
functional formalism using large polarized Gaussian basis sets. The calculated values 
of electronic (agz = 102.8 ) and vibrational (omfe = 4.3 ) polarizabilities indicate 
that the vibrational contribution to the total polarizability of C70 is very small. 

This work is supported in part by the National Science Foundation through CREST 
grant, by the University of Texas at El Paso and by the Office of Naval Research. 
Authors acknowledge the computer time at the UTEP Cray acquired using ONR 
05PR07548-00 grant. 
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